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The spin-orbital model for triply degenerate t2 g electrons on a triangular lattice has been shown 
to be dominated by dimers: the phase diagram contains both strongly resonating, compound spin- 
orbital dimer states and quasi-static, spin-singlet valence-bond (VB) states. To elucidate the nature 
of the true ground state in these different regimes, the model is mapped to a number of quantum 
dimer models (QDMs), each of which has three dimer colors. The generic multicolored QDM, 
illustrated for the two- and three-color cases, possesses a topological color structure, "color vison" 
excitations, and broad regions of resonating VB phases. The specific models are analyzed to gain 
further insight into the likely ground states in the superexchange and direct-exchange limits of the 
electronic Hamiltonian, and suggest a strong tendency towards VB order in all cases. 



PACS numbers: 75.10.Jm, 05.30.-d, 05.50.+q, 71.10.Fd 



I. INTRODUCTION 

Spin-orbital models have become a very active field 
of research in the quest for exotic states of matter and 
novel properties in systems with coupled charge, spin, 
and orbital degrees of freedom. In real materials, each 
class of active orbital, lattice structure, and extent of 
electron filling leads to a different set of problems. An 
overview of this fabulous wealth of possibilities, focusing 
on undoped systems and including a complete review of 
relevant materials, may be found in Ref. 

Among these many models, one of the most exotic 
and mysterious is the triangular-lattice tig spin-orbital 
model considered in Ref. [l|. By this is meant the insu- 
lating system with a single, fully localized electron on 
each lattice site, and occupying one of the threefold- 
degenerate t-2g orbitals, a situation which would be real- 
ized, for example, in undistorted NaTiC>2 or, with holes, 
in C0O2. Depending on which orbitals they occupy, the 
particles have both superexchange and direct-exchange 
interactions whose ratio is a priori unknown, but it was 
showni that different types of bond dimer state of the 
spin and orbital degrees of freedom are favored very 
strongly for all values of this ratio. The resulting, very 
rich, phase diagram of possible states includes candidate 
resonating valence-bond (RVB) states in the limit of pure 
superexchange interactions and, in the direct-exchange 
limit, candidate systems for extremely subtle order-by- 
disordcr selection of special valence-bond crystal (VBC) 
states from very highly degenerate manifolds of dimer 
coverings. 

The question left unanswered despite the extensive en- 
ergetic studies and plausibility arguments in Ref. [l[ was 
whether one could gain more explicit indications for the 
suggested unconventional nature of the true ground state 
in either of these cases. The exotic phenomena which are 
known from simplified models of highly frustrated and 
degenerate systems include true liquid phases with no 
unbroken symmetry in the ground state, topological sec- 



tors and excitations, and fractionalization and deconfine- 
ment of elementary spinon-type quasiparticlcs. Given 
the exhaustive search for both experimental and theo- 
retical realizations of such phenomena in real, electronic 
systems, the triangular-lattice tig spin-orbital model is 
clearly a prime candidate for further investigation. Here 
it can be noted that the preponderance of evidence in 
favor of RVB states suggests further that the candi- 
date (spin-orbital) liquid phase is gapped, and hence 
would have only short-ranged correlation functions and 
massive, fractional spinon-orbiton excitations. The aim 
of this study is to approach the unanswered questions 
through effective quantum dimer models (QDMs). 

A review of QDMs, including their origin, properties, 
and rich associated physics, can be found in Ref. 2|. For 
the present purposes, specifically pursuing the tantaliz- 
ing prospect of an RVB state in a real electronic system, 
the agenda is that laid out as a practical prescription in 
Ref. [3(, namely to map the starting Hamiltonian to the 
triangular-lattice QDM. In the original QDM exposition 
of Rokhsar and Kivelson^ it was shown for the square 
lattice that an exact RVB state exists at one point in 
the phase diagram (henceforth an "RK point"), and it 
was later founds that the same model on the triangular 
lattice allows a rigorous proof for the existence of a true 
RVB phase with gapped spinon excitations over a finite 
regime of parameter space. This therefore offers not only 
a finite chance of any given model having the appropri- 
ate energetics, but also an approach which captures au- 
tomatically the topological criteria associated with RVB 
physics. 

The generic QDM phase diagram^ nevertheless con- 
tains RVB states over at best a rather narrow range of 
parameters. Its bulk is occupied by static "columnar" 
and "staggered" dimer regimes, which respectively max- 
imize or minimize the number of active four-site units, 
or by other possibilities with a more complex breaking of 
translational symmetry on generalized plaquettes. Thus 
the same approach can also be adopted for the direct- 
exchange limit of the starting modelji not in this case 
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with intent to find an RVB phase, but with a view to 
isolating the leading fluctuation term which might be 
responsible for selecting one, or a set of, dimer cover- 
ings which are more energetically favorable, through the 
quantum fluctuations they allow, than all others. 

Because QDMs are in general highly simplified, any 
mapping from a realistic Hamiltonian discards of neces- 
sity many degrees of freedom, and thus there is no sys- 
tematic procedure for their derivation. When proceeding 
from real, S — 1/2, spin-singlet dimers, an immediate 
dichotomy arises between the definition of QDM states 
as mutually orthogonal and the fact that no conventional 
spin-dimcr states are orthogonal (the overlap factor be- 
ing trivially l/y/2 per site). This non-orthogonality 
problem is strongly reduced in the triangular-lattice tig 
spin-orbital model, and entirely absent in its direct- 
exchange limit, for reasons which will become apparent in 
Sec. II. However, this fact also forbids the type of dimer 
overlap expansion exploited in Ref. Q, and mandates 
instead a process which proceeds directly from the elec- 
tronic Hamiltonian. The process of deriving an effective 
QDM from an orbitally degenerate electronic Hamilto- 
nian was followed in Ref. @ , in the course of an analysis 
of LiNiC>2 (a system of e g electrons on a triangular lat- 
tice) , and while the rationale is similar here, the orthogo- 
nality of dimer coverings requires an approach somewhat 
different in detail. 

As a foretaste of the results to follow, in the superex- 
changc limit one finds two different types of "three-color" 
QDM. The individual dimers are either spin singlets with 
three different Savors corresponding to the triplet states 
of the orbitals (ss/ot), or orbital singlets with a conven- 
tional spin-triplet flavor (os/st). In both cases, the 9x9 
t and v matrices of the QDM reflect the strong break- 
ing of translational and rotational symmetry associated 
with the consideration of four-site plaquette units in the 
definition of the QDM. In the former case this drives a 
lifting of degeneracy in favor of the T z = ±1 components 
of the orbital triplet, while in the latter the spin-triplet 
dimers may continue to interchange their flavors. 

In the direct-exchange limit, there is no mixing be- 
tween two pairs of dimers on a plaquette, and both the t 
and v terms of the conventional QDM arc identically zero; 
this rigid locking of dimer color to lattice direction means 
that the effective model reduces to the one-color variant. 
The relevant QDM capturing the lowest-order dimer res- 
onance processes is in fact one defined on six-bond trian- 
gular loops. While this model dictates a specific ground 
state, selected from the extensively degenerate manifold 
of dimer coverings, it does not exclude the possibility of a 
type of one-dimensional physics^ where the character of 
short, fluctuating segments of frustration-decoupled spin 
chains persists despite the high site coordination of the 
triangular lattice. 

Finally, it should be stressed that none of the physics 
discussed here is to be confused with that of the "spin- 
orbital liquid" and "spin-orbital singlet" states intro- 
duced recently in Ref. [7| . The work of these authors de- 



pends intrinsically on a strong spin-orbit coupling inter- 
action \qL-S, and should perhaps therefore be referred to 
more strictly as " (spin-orbit) al" physics. In the present 
study, which follows a line dating back to the seminal 
work of Kugel' and Khomskii^ the terminology "spin- 
orbital" refers only to the connection, dictated precisely 
by the electronic Hamiltonian, between the magnetic ex- 
change interactions and the orbital state of each ion, with 
Ao assumed to be negligible. While this approximation 
is usually taken to be reasonable for 3d ions, Ao may in- 
deed become important when the physics of the system 
depends sensitively on the lifting of high degeneracies. 

The structure of this manuscript is as follows. The 
triangular-lattice ti g spin-orbital model is presented in 
Sec. IL in the form of a minimal review of the contents of 
Ref. [l[ required as a basis for the current analysis. The 
orthogonality of electronic states with different spin or 
orbital colors (colored dimers) is demonstrated and the 
derivation of appropriate QDMs outlined. In Sec. Ill, the 
generic multicolored QDM is defined and its properties 
are studied, illustrating the circumstances under which 
the color degree of freedom leads to additional topologi- 
cal sectors, possibilities for RK points and RVB ground 
states, and "color vison" excitations associated with the 
topological properties. Section IV returns to a detailed 
analysis of the specific QDMs deduced for the three cases 
of most interest here, namely (ss/ot) and (os/st) dimers 
in the superexchange limit and spin-singlet dimers in the 
direct-exchange limit. A short discussion and summary 
are presented in Sec. V. 



II. MODEL 

The derivation begins with the Hubbard Hamiltonian 
for tig electrons hopping on a triangular lattice, usually 
realized on the (111) planes of a cubic, perovskite-based 
structure whose local symmetry preserves the orbital de- 
generacy The geometry of the system is summarized 
in Fig. 1, and the average filling of 1 electron per site 
corresponds to Ti 3+ or V 4+ ions. To avoid undue repeti- 
tion of published material, only the minimal numbers of 
equations and explanations required for a self-consistent 
presentation are included here; the reader wishing a com- 
plete exposition is requested to consult Ref. [l| ■ 

Specializing immediately to a magnetic insulator with 
on-site repulsion U much greater than the diagonal and 
off-diagonal (in orbital "color") electronic hopping inte- 
grals, respectively t' e and t e , the site occupation is strictly 
m a + riib + rii c = 1. The orbitals have a direct connec- 
tion to the specific lattice directions, and in Fig. 1 the 
color red for the d yz orbital, is associated with a, green 
(d X z) with 6, and blue (d xy ) with c. At second order in 
a perturbative treatment, the magnetic Hamiltonian can 
be written in the schematic form 

H = J S T~L S + J m T~L m + JdUd, (1) 
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FIG. 1: (color online) (a) Schematic representation of hop- 
ping processes for t2 g electrons in triangular geometry which 
contribute to magnetic interactions on a representative (c- 
axis) bond (ij). The orbitals are represented by different col- 
ors (greyscale intensities). Superexchange processes involve 
O 2p z orbitals (violet), and couple pairs of a and b orbitals 
(red, green) with effective hopping elements t e , interchanging 
their orbital color. Direct exchange couples c orbitals (blue) 
with hopping strength t e . (b) Pairs of ti g orbitals active in su- 
perexchange and (c) single orbitals active in direct exchange; 
horizontal bonds correspond to the situation depicted in panel 
(a). 
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are dictated by the different possible energy states of the 
virtual d\d^ ^ d^d® excitation, which are determined by 
the Hund coupling Jh = f]U . In real 3d transition metal 



ions, r\ is of order 0.10-0.15. Under a local transformation 
which takes account of the color-exchanging nature of the 
superexchange hopping processes, 
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where the scalar product in Ay is the conventional ex- 
pression for pseudospin-1/2 variables and the cross prod- 
uct in B.jj is defined as 
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Contributions of the form B, 
(U). In these expressions, n v 
superexchange-active electrons (those able to hop on the 
bond in question) and the number of superexchange- 
inactive electrons (which, however, are active in direct- 
exchange processes): thus the pseudospin scalar and vec- 
tor product interactions are relevant only when the elec- 
trons on both sites are active, while situations with only 
one active orbital still give nontrivial, spin-dependent 
contributions. 

It is clear that the superexchange contribution on a sin- 
gle bond is minimized either by an orbital-singlet, spin- 
triplet state, or by a spin-singlet, orbital-triplet state 
[henceforth (os/st) and (ss/ot)]. Their energies 
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are degenerate for r\ = 0, while the (os/st) state is favored 
for finite rj. Because the superexchange hopping term is 
off-diagonal, the orbital singlet is the state 
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in the original electronic basis, while the orbital triplets 
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Contributions from direct-exchange processes take the 
form 
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where it is clear once again that far the most favorable 
energies are obtained from dimer states, but only those 
creating spin singlets from two electrons with the bond 
color. Expressions for H m can be found in Ref. [l|, but, as 
explained there, the very different nature of the two types 
of singlet state mean that H m has a qualitative effect on 
the physics of the system only in rare situations. 

Before analyzing the model in more detail, it is neces- 
sary to delimit the parameter space of interest to this 
study, and also to make one essential and completely 
general comment concerning dimer states of colored elec- 
trons. First, the focus of this manuscript is only the two 
limits of pure superexchange and purely direct exchange. 
In the superexchange limit, only the regimes with Hund 
coupling ratio 77 = and with physical values rj ~ 0.1 
(which is "large" in a sense to be made explicit below), 
will be considered. In the direct-exchange limit, only 
i] = is of interest in the context of QDMs: the physics 
of dimer-bascd models is determined by purely static 
contributions at any finite 77^ and changes away from 
dimer-based states only at the unrealistically large value 
i] = 0.2. Thus there are three situations to elucidate, 
all of which have already been shown to be dominated 
by dimer formation. This restriction is made partly for 
reasons of practicality, partly to concentrate on the essen- 
tially different physics of the differing limits, and partly 
as a consequence of the above observation concerning the 
irrelevance of H m and absence of meaningful mixed pro- 
cesses. 



A. Effective dimer overlap 

As noted above, a serious problem in mapping real 
spin models to QDMs is the "non-orthogonality catas- 
trophe," the fact that all dimer coverings |c) have finite 
overlap with all others. This is a significant impedi- 
ment to many forms of both analytical and numerical 
progress. For S = 1/2 entities forming bond spin-singlet 
states, the generic overlap of two singlets at the same 
site is simply the prefactor of the singlet wave function, 
a = l/v2- In the present context, the most mean- 
ingful calculation involving a is to consider the over- 
lap of the wave functions of two "horizontal" and two 
"vertical" SU(2) spin singlets on a four-site plaquette, 
where it is easy to show from the common elements of 
\h) = a 2 (|lt'2|> - |1| 2t»(|4t 3|) - |4| 3t» and 
\v) = a 2 (|4t U) - \U lt»(|3t 2|) - |3| 2t» that 
o = (h\v) = 2a* = 1/2. 

The single most important generic feature of raising 
the number of dimer colors, noted already in the con- 
text of 1/N expansions^ is expected to be that a is re- 
duced, perhaps by a power. However, this hypothesis is 
dependent on the model, and for a model such as the cur- 
rent one, where electron orbital color sectors change with 
bond direction, the overlap matrix is not uniform, and in 
some cases is identically zero. For the three cases of in- 
terest here, this can be shown very simply by repeating 



the exercise above for (os/st), (ss/ot), and (ss/77) units, 
where 7 indicates the bond color. For each of the two 
spin-orbital triplet cases there are nine possible overlap 
matrix elements to consider, not all of which will be il- 
lustrated explicitly. Taking a plaquette of parallel a and 
c sides, one has 



h++ (os/st)) 


= a 2 (|l?-t 2rf) - 


|l<7t 2fft)) 




x(|4rt 3rf) 


- |4flft 35t)), 


\v++ (os/st)) 


= a 2 {\Ag\lg\)- 


|4&tl6t)) 




x(|3flt 2fit) 


-W 2H», 


|/ioo(os/st)) 


= a 4 (|lrt2r|) + 


114 24) 




-llflt 2.4) 


-\ig±2fft)) 




x(|4rt 34) 


+ 144 34) 




-|4flt3<?4) 


-|4ff|3gt), 


1 woo (os/st)) 


= 4(|4.4 1.4) + \4gl Ifft) 




-mm) 


- m 




x(|3flft 2.4) 


+ |3<?+ 2gt) 




— |36t 260.) 


- |3fcd. 2&t», 


h++ (ss/ot)) 


= a 2 (|M2.g|)- 


|142flt)) 




x(|4rt Sgi) 


-144 3.4)), 


[«++ (ss/ot)) 


= a 2 (|4gtl6|)- 


14.414)) 




x(|3flt 26|) 


-13.424)), 


|/ioo(ss/ot)) 


= a 4 (|lrt2r|) - 


11424) 




+ list 2s;) 


-1 1.4 24)) 




x(|4rf 34) 


- 144 34) 




+|4flt 3<?i) 


-14434)), 


1 woo (ss/ot)) 


= a 4 (|4gt lgl) - 


K4 1.4) 




+|4&t 14) 


- 144 14)) 




x(|3fft 2<?4) 


- 13.4 24) 




+|36t 2&J.) 


-13424)), 


/i(ss/cc)) 


= a 2 (|16t2fe|) - 


I1424) 




x(|46f3fei) - 


-144 34)), 


u(ss/aa)) 


= a 2 (|4rt 14) - 


44 14) 




x(|34 24) 


-134 24)), 



and similarly, whence 
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o(ss/ 7 7) = 0. (15) 
For (os/st) dimers there are finite overlap matrix ele- 
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FIG. 2: (color online) QDM representations, (a) Three-color 
QDM as obtained in the superexchange limit for both (ss/ot) 
and (os/st) dimers; dimers of any color may occupy any bond, 
but not all processes (t or v) may be allowed on plaquettes 
containing two dimers. Dimer colors do not correspond to the 
colors in Fig. 1. (b) QDM obtained in the direct-exchange 
limit for (ss/77) dimers, whose colors do correspond to those 
in Fig. 1. Bonds of the three different orientations have only 
one color each, and thus the color is no longer a degree of 
freedom. This defines an effective one-color QDM. 

merits within each S z sector, all (but one) of considerably 
higher order than the conventional spin-singlet problem. 
However, the 9x9 overlap matrix is rather sparse. For 
(ss/ot) dimers, the orbital color combinations forbid any 
overlap matrix elements other than for pairs of T z = 
dimers. For bond-colored, spin-singlet dimers there are 
no overlap matrix elements at all. Thus the highly non- 
uniform nature of the matrix elements in the superex- 
change model and the fact that they vanish in the direct- 
exchange model preclude any sort of meaningful expan- 
sion in the inverse of the overlap matrix^ As noted in 
Sec. I, the matrix elements t and v of the relevent QDMs 
will instead be deduced directly at second-order in per- 
turbation theory using the electronic Hamilitonian. 



three analogous states for a different lattice direction are 
composed of a different t2g~color pair. 

However, for physical values of rj around 0.1, the 
(os/st) dimers are preferred as the direct lifting of de- 
generacy in the one-bond Hamiltonian outweighs the ki- 
netic effects. This regime is delimited by a lower bound, 
which for the four-site cluster is rji — 0.03, and an upper 
bound of i] u ~ 0.15 where the ground state becomes fer- 
romagnetic by a lifting of the spin-triplet degeneracy,! 
Within this range, one has a three-color QDM similar to 
the (ss/ot) case [Fig. 2(a)], but one in which the colors 
correspond to the spin triplet states of the two electrons, 
S z = 1, 0, and —1. The triplet states in the spin sector 
are the same for every bond direction, and this fact is 
responsible for the differences from the (ss/ot) situation 
observed already in the overlap matrix elements. The 
QDM matrix elements can also be expected to differ in 
several respects. 

Summarizing the situation to this point, in the su- 
perexchange limit of the model one has a number of ener- 
getic criteria which can be compared to the RVB criteria 
of Ref. [||. The former arei (i) a very strong tendency 
to dimer formation, (ii) a large semi-classical degeneracy 
of basis states, namely the set of coverings formed from 
these dimers, and (Hi) that resonance processes on four- 
site plaquctte units provide a significant energetic con- 
tribution. From (Hi), a minimal triangular-lattice QDM 
should in principle contain, within the dimer-flipping 
term t and the static pair term v, the leading correc- 
tions to the static VB energy (determined by J s ). The 
question is then reduced to whether the regime t/v > 1 
of RVB physics^ can be attained. 

C. Direct— exchange model 



B. Superexchange model 

In the limit 77 = 0, the bond Hamiltonian of Eq. ^ 
takes the more transparent form 

+ \{^i + n r<)-\- (16) 

The six-fold single-bond degeneracy of (ss/ot) and 
(os/st) dimers is lifted in favor of the former by quantum 
fluctuations on the off-dimer bonds at r\ = 0^ leaving a 
threefold degeneracy corresponding to the orbital triplet 
states of the two electrons, T z = 1, 0, and —1. Thus 
one expects a QDM with three different flavors of dimer, 
a situation depicted in Fig. 2(a). Note that these "col- 
ors" are not the same as the spectral colors in Fig. 1, 
which are associated directly with the symmetry of the 
t2 g levels, but correspond to three possible orbital states 
on each bond of a given direction [Eqs. (fTTj) - (fl"3"|) ]. The 



Here only the r\ = limit of the model, 
Udo = Jd^2 ~ jj n ii n 3i ( 17 ) 

('0')ll7 

-4(^7(1 - n ii) + (1 - n n ) nji y 

will be considered. As noted above, the spin-singlet state 
optimizing the bond energy can be formed only from two 
electrons of the bond color. Thus the three apparent t^g 
colors of the dimers are locked rigidly to the bond direc- 
tion, with no finite matrix elements for directional fluc- 
tuations, a situation depicted in Fig. 2(b). The effective 
QDM will then be equivalent to a one-color model. 

In addition to the studies of Ref. [l[ , this model has also 
been considered in some detail by Jackeli and Ivanov^ 
These authors demonstrated that the extensively degen- 
erate manifold of (bond-colored) dimer coverings forms a 
set of exact eigenstates of the second-order Hamiltonian. 
In deducing the leading perturbations which would select 
a type of order, they allowed finite values of Jh, which 
lead at order rf to a ferromagnetic interaction on inter- 
dimer bonds with one (but not two) dimers of the bond 




FIG. 3: (color online) Representation of allowed dimer fluc- 
tuation processes in the direct-exchange limit, shown for a 
system of only one dimer color. These processes define loops 
whose sides have even numbers of bonds, and which (to main- 
tain complete dimer covering of the lattice) must enclose an 
even number of sites. Illustrated are the lowest three loops, 
which have lengths of 6, 10, and 14 bonds. 



color. As a result, particular static dimer configurations 
are selected (the type of state is unique, albeit with a 
translational and rotational degeneracy of 60), but there 
remains no contribution from the positional resonance of 
dimers. The aim here is to investigate the consequences 
of the lowest-order dimer resonance processes of a type 
which would give valid contributions to a QDM. 

In this limit of the model, electrons in dimer sin- 
glets may hop only in the direction corresponding to 
their bond color. Thus virtual processes are highly re- 
stricted and the true dimensionality of the system is far 
lower than the connectivity of the triangular lattice would 
suggest.— Because a four-site plaquette contains no two 
bonds of the same direction, not only is the overlap ma- 
trix element zero (above) but the t and v terms of the 
conventional, triangular-lattice QDM are also identically 
zero. A QDM based on higher-order terms will therefore 
be required. The only kinetic contributions on the trian- 
gular lattice are given by paths of even side-lengths, and 
which enclose even numbers of sites. Some valid loops 
are shown in Fig. 3. Here the contributions t' and v' 
from fluctuation processes on six-bond triangles will be 
computed, while the leading perturbations from higher 
allowed loops (the 10-bond trapezium and 14-bond ir- 
regular pentagon) can be argued to be small. The ground 
state of the resulting t!-v' QDM should contain the lead- 
ing effects of dimer resonance processes on the highly de- 
generate manifold of static VB coverings^ One may then 
also consider whether the energy of this state is more fa- 
vorable than the spontaneous one-dimcnsionalization of 
the system into effectively decoupled Heisenberg chains 
of a single orbital color. This outcome, the lowest-energy 
possibility found in Ref. [l| , is generally held to be rather 
unlikely, although from the ID nature of the electronic 
processes in this limit of the model it is not entirely im- 
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FIG. 4: (color online) (a) Dimer covering of the two-color 
QDM on the triangular lattice, (b) Representation of the 
eight possible configurations of two colored dimers on a flip- 
pable (ac) plaquette. 



plausible. 

Commenting briefly on further models which may arise 
within the t2 g system on the triangular lattice, clearly 
there are many of these even at unit filling per site. The 
entire range < J' / J < oo is of obvious physical im- 
portance, and in principle presents a model with four 
types of dimer. However, the lack of symmetry between 
the two limits implies automatic selection effects which 
will favor a smaller number of these dimers (those max- 
imizing virtual kinetic proceses) . The results of Ref. [l[ 
suggest a first-order transition from the superexchange 
regime (QDMs with up to three colors) to the direct- 
exchange regime (a one-color QDM) at an intermediate 
value of J' I J but with no intermediate phase. A sec- 
ond possibility is the regime of very high rj, where for all 
values of J' /J the preference for dimers will melt, and 
one obtains very highly degenerate sets of color configu- 
rations among which the (now subdominant) Heisenberg 
interactions are responsible for selection effects. 

To conclude this section, specific colored QDMs can 
be obtained for the three cases of interest and, within 
the approximations inherent in the derivation process (by 
which is meant the breaking of translational symmetry 
and the discarding of very many degrees of freedom) , all 
are expected to yield additional information concerning 
the physics of the ground state of the starting electronic 
Hamiltonian. Before proceeding to a detailed analysis 
of these models (Sec. IV), it is necessary to consider the 
physical possibilities contained within QDMs possessing 
an additional color degree of freedom. 
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III. GENERAL COLORED QDM 

Consider a QDM on the triangular lattice with two 
different types of dimer, which will be labeled by the 
colors red and blue. An arbitrarily chosen dimer covering 
is shown in Fig. 4(a). To systematize the description of 
this system, one begins by noting that there are eight 
possible states on each rhombic plaquette: taking an (ac) 
plaquette for illustration [Fig. 4(b)], let these be denoted 
by i = rr, rb, br, and bb for pairs of a-axis dimers and 
(rr), (rb), (br), and (bb) for pairs of c-axis dimers. The 
QDM Hamiltonian has the form 

3 2N 

# = EE^ ( 18 ) 

7—1 2—1 

where there are 2N plaquettes of any one orientation for 
a system of 2N sites (N dimers), and 7 labels the three 
plaquette orientations on the triangular lattice. The pla- 
quette Hamiltonian is given by 

H p = -trr,rr(\& , )(&\ + \&){&'\) 

-Ur,rb{\&){&\ + \ff)(&\)+... 

-UbM{\&)(&\ + \&)(&\)+--- 

+V„,rr(\JS')( i ff\ + \&)(&\) 
+Vrr,r b (\£')(£'\ + \&)(&\) +... 
+V b r,bb(]&)(&\ + \&){&\)+--- 

and thus is represented by an 8x8 matrix acting on the 
state vector [rr, rb, br, bb, (rr), (rb), (br), (bb)]. This ma- 
trix has a "diagonal block" and off-diagonal blocks Uj 
and t*y The 4x4 matrix specifies the static energies 
of the rr, rb, br, and bb configurations, and the assump- 
tion of plaquette symmetry between a- and c-axis dimer 
pairs, enforced here, means that this symmetric matrix 
also defines the analogous energies for (rr), (rb), (br), 
and (bb). The 4x4 iy matrix specifies the flipping ener- 
gies of rr, rb, br, and bb into (rr), (rb), (br), and (bb), and 
conversely; the symmetry = t(ju (Hcrmitian conju- 
gation) ensures that is a symmetric matrix and the 
symmetry = tajj (one-plaquette axis permutation) 
ensures further that tij is symmetric. Thus the brackets 
are redundant in the labels of Vij and tij . 

The nature of the two-color triangular-lattice QDM 
depends on the matrix elements of Vij and tij, which 
could in principle vary between the limits 

= v[l, 0,0,0; 0,0, 0,0; 0,0, 0,0; 0,0, 0,1], 
Uj = t[l, 0,0,0; 0,0, 0,0; 0,0, 0,0; 0,0, 0,1], (19) 



and 

Vij = v[l, 1,1,1; 1,1, 1,1; 1,1, 1,1; 1,1, 1,1], 

hj = t[l, 1,1,1; 1,1, 1,1; 1,1, 1,1; 1,1, 1,1]. (20) 

The system represented by Eqs. (|19j) is one where po- 
tential and kinetic energy are gained only on plaquettes 
containing two dimers of the same color, and hence this 
is a dilute version of the one-color QDM. While some ad- 
ditional symmetries may arise (below) from the conser- 
vation of each color individually, the energetics and dy- 
namics of the system appear to be somewhat trivial as a 
result of the dilution. Equations (|20|) represent a "color- 
blind" system where the standard QDM processes occur 
with the same amplitudes irrespective of dimer color, and 
is therefore equivalent to a one-color QDM but with a 
redundant degree of freedom. While it is the signs of 
the elements of and t^ which are most important in 
determining the ground state of the model, other than 
Vu there is in general no need for these to be real, and 
one could consider matrix elements which have different 
signs or phase factors occurring in different patterns. 



A. Topological structure 

One of the most important properties of QDMs is the 
fact that they display a topological structure (also "topo- 
logical order" )j£ which is characterized by the presence of 
topological invariants and associated gauge symmetries. 
For a colored QDM it is therefore natural to investigate 
whether there exists some additional, nontrivial topologi- 
cal structure arising from the presence of the color degree 
of freedom. 

For this one must first recall the origin of the Z2 gauge 
symmetry which exists in the conventional QDM. This is 
ultimately a consequence of the condition that each site 
must contain only one hard dimer on one of the bonds 
connected to it. If one considers the number of dimers 
intersecting an arbitrary line which spans the system, it 
is this constraint which is responsible for the fact that 
any dimer-rearrangement process (which defines a loop) 
can change the number of dimers intersecting the line 
only in multiples of 2. This property cannot be affected 
by any set of color matrix elements tij and within 
the framework presented above, and therefore this Z2 
symmetry cannot be embedded within some higher sym- 
metry group. Alternatively stated, while it is possible 
to alter dimer colors without altering the number which 
intersect the arbitrarily chosen line, it is not possible to 
alter the number intersecting the line without altering 
the numbers of dimers of each color which intersect it. 

Thus the symmetry of the colored QDM can be at best 
Zi®G, where G is a different symmetry group. A number 
of possibilities could be considered which avoid this hard 
constraint, including "soft" dimers which must only have 
unit squared amplitude on the sum of the bonds at a site, 
multiple dimers per site (for example a representation of 
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a S = 1 spin system), a diluted model with undimerized 
sites or pairs of sites, and models of quantum trimers or 
quantum quadrumers on particular lattices and with par- 
ticular matrix elements analogous to Ay and Wy. While 
some of these may indeed be constructed to deliver topo- 
logical structures and gauge symmetries beyond Z2, they 
will not be considered further in the present context. 

Returning to a consideration of the number of dimers 
intersecting an arbitrary line across the system, the quan- 
tity ./V = TV,, + Nf, (the total numbers of red and blue 
dimers in the system) is always a constant, while n r + 
(the numbers of red and blue dimers intersecting the line) 
may change only in units of 2. This property defines the 
four topological sectors on a torus and is reflected in the 
Z2 symmetry discussed above. 

Nontrivial color sectors are characterized by n r — rib- 
This quantity can vary from — TV to +TV in steps of 2; 
there are TV sectors if TV is odd and TV + 1 if TV is even. 
Thus a trivial color-sector structure, where sectors with 
(N r ,Nb) = (K,L) are connected to sectors with (K ± 
1, L=pl), is characterized by (n r — rib)/2 changing in steps 
of 1. Two nontrivial alternatives arise in the two-color 
QDM. 

(i) If the matrices Ay and i>y have the forms 
[a, 0, 0, 0; 0, b, c, 0; 0, c, b, 0; 0, 0, 0, d], TV,, and TV fc are 
conserved individually and {n r —rib)/2 is a constant. 
Thus there are TV or TV + 1 topological color sectors 
which cannot be mixed. The symmetry may be la- 
beled Z2<E>Z. Models with an infinite set of discrete 
topological invariants labeled by the set of integers 
Z constitute a realization of an ensemble of string 
nets, the theory of which can be found in Ref. [Ioj |. 

(ii) If the matrices Ay and i>y have the forms 
[a, 0, 0, e; 0, b, c, 0; 0, c, b, 0; e, 0, 0, d\, terms intercon- 
verting between two red and two blue dimers are 
also permitted. In this case, sectors (K, L) are con- 
nected to sectors (K ± 2,L =p 2) and (n r — rib)/ 2 
(mod 2) is a topological invariant. Sectors (K, L) 
and (K±l, L^pl) have no overlap, and changing the 
color of a single dimer is a process which results in 
a topological defect (below), one which cannot be 
rectified without a string of local processes which 
extends to the boundary of the system. In this case 
one has four additional color sectors on a torus, and 
the symmetry is Zi®Zi'. one Zi is for the hard 
dimer structure and the other for the color struc- 
ture. 

If further elements in case (ii) were to be made finite, 
then the Hamiltonian matrix would mix sectors differing 
by (n r — rib)/ 2 = 1, and the trivial topological structure 
would be restored. 

Thus the presence of a dimer color in a QDM can lead 
to additional topological structure in the color sector. 



B. RK points and RVB phases 

Another key property of QDMs is the fact that they 
display Rokhsar-Kivelson (RK) points. A further essen- 
tial question for colored QDMs therefore concerns the 
existence of RK points in models described by the differ- 
ent types of matrix above. The RK point of the original 
QDM occurs at A = v. It has the property that the 
equal-amplitude superposition (with equal phase) of all 
possible dimer coverings in each topological sector, a "liq- 
uid" state with no local order, is the ground state. While 
on the square lattice there are gapless excitations at the 
RK point, on the triangular lattice this ground state is 
completely gapped and the properties of the RK point 
are preserved across a phase of finite extent. This is the 
RVB phase. 

The equal-amplitude superposition \ip) = -j= Y)\ " |c), 
where |c) represents the TV C dimer coverings in the cho- 
sen topological sector, has energy (ip\H\ip) = (v — t)(n ), 
where n fl is the number of flippable plaquettes in a cov- 
ering. To show that this wave function is an eigenstate 
requires writing the QDM Hamiltonian as a sum of pla- 
quette projectors, which in turn requires the condition 
V = t. At the RK point, the energy is zero, while on 
the triangular lattice it is negative over a range of val- 
ues 1 < t/v < 1.3, beyond which there are transitions 
to different ground states (of ordered dimers or groups 
of plaquettes, hence breaking the translational and rota- 
tional symmetry of the lattice). 

The proof that a more general \ip) is an eigenstate of a 
more general QDM Hamiltonian sets the restrictive con- 
dition that H must still be expressible as a sum of pro- 
jectors. In a two-color model of the type considered in 
this section, where H contains many terms Ay and Uy, 
the RK condition is that Ay = uy for each matrix com- 
ponent. The wave function in this case contains all cov- 
erings, in each topological sector (where now the color 
may increase the number of sectors), of dimers of both 
colors. If this structure is not maintained, \ip) will not 
be an eigenstate, but may still be the ground state, and 
will retain a gap to all excitations up to the boundary of 
the RVB phase. 

The condition that \ip) be the ground state of H is, 
up to a point, less strict. The energy of the two-color 
version of the wave function \ip) is 

E = ^9ij(vij -Uj), (21) 

where gy is a statistical factor accounting for the proba- 
bility of finding each specific plaquette state (here those 
plaquettes not only containing two dimers but also with 
a finite matrix element in Ay and/or uy for processes 
involving these). In fact <?y = (nfj) is the average num- 
ber of flippable plaquettes of each type in the equal- 
amplitude superposition. For equal dimer numbers, 
TV,. = TVb, in a two-color model, all coefficients <?y are 
equal, but this is no longer the case if, for example, 
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in the \/l2 x VT2 formation^iii This phase-separation 
effect becomes less pronounced as the values of the other 
matrix elements in Hij (for example contributions from 
flippable red-blue plaquettes) are increased. 




FIG. 5: (color online) Schematic representation of a two-color 
dimer system undergoing phase separation: if t rr ,rr > v rr ,rr 
but tbb,bb < Vbb.bb, a dynamically driven phase separation will 
occur (mitigated by larger values of matrix elements mixing 
with local rb and (rb) states). The "staggered" configuration 
of the predominantly blue region is static, whereas the red 
region undergoes local, symmetry-restoring resonance pro- 
cesses. 



N r =/= Nf,, or if a more complex model defined on a more 
complex geometry is considered. 

The physics of the colored QDM is in this regard 
not qualitatively different from that of the conventional 
QDM. The equal-amplitude superposition does not gain 
a large amount of energy on each plaquettc, but it is 
the ground state because it profits from quantum fluctu- 
ations on every single plaquette, and not merely from a 
symmetry-broken subset of these (which can be at most 
1/6 of the total plaquette number). The freedom to have 
every component obey 1 < Uj/vij < 1.3 clearly defines 
a relatively broad regime of parameter space: should the 
ratio stray outside these bounds for some components, 
the preference of the remaining components for full res- 
onance will act to maintain the RVB phase. 

However, in this situation a multicolored QDM is sus- 
ceptible to phase separation. As a specific but simple 
example, consider a model with only the terms v rr ^ rr , 
Vbb,bb, trr,m and tbb.bb (all positive) large, while other 
matrix elements arc finite [to avoid the dilution paralysis 
problem of Eq. (|T9")) ] but very small. The large elements 
obey t rr ,rr > v rrtrr and tbbbb < Vbb.bb, whence even if 



> VbbM — tbb.bb, 



(22) 



so that the sum of kinetic terms exceeds that of poten- 
tial terms, flippable blue plaquettes will entail a local 
energy cost. This can be avoided if blue dimers tend to 
adopt a "staggered" phased meaning any one of the set 
of dimer coverings ensuring that there are no flippable 
blue plaquettes. The transition to the staggered phase in 
the conventional triangular-lattice QDM is of first order, 
making the tendency towards phase separation strong. 
In the meantime, red plaquettes may continue to profit 
from resonance across a fraction of the system consistent 
with the ratio N r /N (Fig. 5), but the wave function is 
now far from an equal-amplitude superposition of all pos- 
sibilities. At the opposite boundary of the RVB region, 
one of the colors could be expected to begin to cluster 



C. Color visons and colored Majorana fermions 

The presence of topologically distinct Z2 color sectors 
in case (ii) above results in color vison excitations. While 
local processes contained in the matrix elements of tij 
and Vij result in excited states in the same topological 
sector, some states with very minimal differences can be 
in different topological sectors. An example in case (ii) 
is the exchange of color on a single dimer, a process not 
permitted by the local matrix elements. The (low-lying) 
excitation which results is topological in nature: to con- 
nect the initial and final states requires a string of local 
processes which spans the system^ An alternative, sim- 
ilar to an intermediate stage in the process of "repairing" 
the topological defect, is the presence of defect pairs, or 
color visons connected by a finite-length path. These de- 
fects are exactly analogous to the conventional Z2 visons 
which arise from the dimer positions, but here they are 
(for a pure color vison, and not a vison of mixed charac- 
ter) present only in the color sector, i.e. without altering 
any dimer positions. Topological excitations will also be 
present in the Z color sector, for which the reader is re- 
ferred to Ref. [13]. 

Statements about color visons may be made rigorous 
by formulating a Majorana-fermion representation of the 
statistical average over the different dimer coverings. Fol- 
lowing Ref. [l2T ] , the complete ensemble of dimer configu- 
rations may be represented by the partition function of a 
set of real, auxiliary fermionic variables a™ correspond- 
ing to each lattice site, 



Z = Jl[ da\ n cxp[ Y, afA^r'af] 



Pfaff(^ m ) 



(23) 

where the right-hand side is the Pfaffian of the matrix of 
effective Majorana-fermion hopping amplitudes on each 
bond (IV) of the lattice. In the two-color QDM, two 
colors of Majorana fermion arc required for each site, and 
this is represented in the index m. The matrix elements 
A™, m take the values ±1 for neighboring sites I and V, 
and otherwise. Hence each bond is described by a 2x2 
matrix 



Arr 
IV 
Abr 



Arb 
A W 
Abb 



(24) 



which is net antisymmetric under exchange of I and V 
(A™, m = -A™ m , A™, m ' = -A™[ m ). Because all but the 
sparsest models contain plaquette processes which on any 
given bond can exchange red and blue dimer ends, and 
change the number of dimer ends of each color, all matrix 
elements are finite in the (r, b) basis. The condition which 
ensures that all dimer coverings appear in the partition 
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function with the correct relative sign is that the product 
J| r A u> = — I on any loop Y consisting of an even 
number of bonds {IV): the product of an even number 
of 2x2 matrices around a plaquettc (of 4 to 2N sites) 
should be the negative 2x2 identity. The two —1 terms 
in — 7 correspond to the two Z2 sectors in the two-color 
QDM. 

The Pfaffian formulation becomes more transparent on 
symmetrizing the Majorana fermions, 

of = ^K + af), 

of = ^(o[-aJ), (25) 

whence A\f, = (A$f,+A$+A§+A$)/2 and similarly for 
the other three elements of Aw in the symmetrized basis. 
The situation is particularly simple in a color-symmetric 
problem, by which is meant one with N r = Nb, v rr ^ rr = 
Vbb.bb, and t rTirr = tbb,bb] other situations are qualitatively 
similar but notationally more complex. In the color- 
symmetric case, the model is invariant under interchange 
of the two colors, r O b. Hence A\[, = Aj$ = A™ m and 

A W = A w = A w m ' > which leads to 

(Aft, A&\ = (AIT+Affr' \ 

{Aft Afp) \ o Ayr -Ayr')- 1 ' 

Thus the bond matrix is diagonal in the {r+b, r—b) basis. 
The plaquette product then reduces to two products of 
even numbers of scalars, each of which has value — 1. The 
upper product corresponds to the (hard-dimer-relatcd) 
Z2 symmetry of the conventional one-color model and 
the lower to the Z2 symmetry of the color sector. 

Defective plaquettes, those where the product is no 
longer —1, describe states which do not belong in the 
same topological sector, and hence are topological de- 
fects. Defects in the upper index of the symmetrized 
bond matrix are conventional visons, while those in the 
lower index are color visons. While it is difficult (for ob- 
vious reasons) to represent a single vison in a finite region 
of a system, Fig. 6 presents a schematic representation 
of two color visons connected by a path of finite length. 
The Z2 symmetry of the color sector dictates that a vison 
is its own antivison. 

The discussion can be extended in a straighforward 
manner to three dimer colors. It is important to note 
first that three colors do not imply a Z% symmetry, or 
indeed anything similar to this. A Z3 symmetry would 
require three separate topological sectors for each peri- 
odic direction on a cylinder or torus, and no aspect of 
the purely pairwise interactions in a QDM permits terms 
creating such a structure; this is not a question of the 
choice of the matrix elements in t ij and Vij , but one con- 
cerning the architecture of the model. With green (g) as 
the third color, the analogs of cases (i) and (ii) above are 
rather the following. 

(is) For interactions conserving the number of dimcrs 
of each color, there are N(N + l)/2 (for odd N) 




FIG. 6: (color online) (a) Arbitrarily chosen dimer configura- 
tion in a two-color QDM. (b) A color vison is introduced in 
the lower left corner by exchanging the color of a single dimer 
(here from red to blue) without altering the positions of any 
dimers of either color. Because this configuration would not 
be in the same topological color sector as that in panel (a), 
attempting to restore configuration (a) by applying the lo- 
cal matrix elements of case (ii) results in a path which spans 
the system. If, however, a second color vison is also present, 
in the upper right corner, then the vison pair may annihi- 
late. For simplicity, the visons are represented as residing on 
flippable plaquettes in which the color-exchanged dimer is in- 
volved; dealing with color-exchanged dimers not initially on 
a flippable plaquette is a small extension, and more general 
considerations can be used to isolate visons as triangle-based 
entities.— (c) One possible path of local processes connecting 
the two color visons. 

or {N + 1){N + 2)/2 (for even N) separate and 
unmixed color sectors which would be labeled by 
Z®Z. 

(I13) For interactions conserving the number of dimer 
colors in each pair modulo 2, there is a set of 
"three" separate Z2 subsectors corresponding to 
(n r — nb)/2, (rib ~ n g)/2 and (n g — n r )/2 (mod 
2). One of the three constraints is trivially redun- 
dant (a necessary consequence of the other two), 
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and hence the model has a Z^®Z2®Z2 topologi- 
cal structure with one hard-dimer sector and two 
independent Z2 color sectors. 

A Pfaffian formulation of the statistical average over all 
dimer coverings would require three Majorana fcrmions 
per site. In case (ii3 ) it would contain one conventional 
vison plus two types of rgb vison. 



IV. TRIANGULAR-LATTICE t 2g MODELS 

In this section, let J = Atl/U denote the bare superex- 
change bond strength (J s in Sec. II) and J' = At'^/U the 
bare direct-exchange interaction ( J c i in Sec. II). Some of 
the results to follow lie beyond the mean-field Ansatze 
used in Ref. [l|. All derivations will be restricted to the 
lowest relevant order in the electronic Hamiltonian. Fur- 
ther matrix elements would naturally be present at higher 
orders, and although these are small, they may have qual- 
itative effects such as the mixing of color sectors (i.e. the 
introduction of further off-diagonal elements in the QDM 
matrices). 

In the results of this section, the QDM matrix elements 
v and v' are always negative because all of their even- 
order contributions are net energy gains. The signs of t 
and t' depend on the choice of phase for each singlet, and 
hence on an arrow-direction convention applied to the 
triangular lattice. These signs are, however, unimportant 
in the analysis of a QDM, where they are a matter of 
convention (easily altered by a change of bond phases) 
and thus will not be tracked explicitly. In fact the signs of 
t and if here, with no further manipulation, are generally 
positive (bearing in mind the overall minus sign in the 
QDM definition) , because all loops involve even numbers 
of plaquette edge bonds of the same type. The sign of the 
ratio t/v is, however, crucial in determining the ground 
state. 



A. Superexchange limit, r\ — 

Here the ground state of the system is composed of 
(ss/ot) color triplets. As noted in calculating the over- 
lap of different dimer states in Sec. II, this model has 
quite asymmetric matrix elements due to the breaking 
of translational symmetry involved in considering four- 
site plaquette units. This is due specifically to the direct 
constraint on the orbital space by the choice of plaquette, 
and hence is stronger for (ss/ot) than for (os/st) dimers. 

Consider an (ac) plaquette, as represented in Fig. 7(a). 
Let the label + denote the state \rg), — the state \gr) 
and the state fyrr) + \gg) [see Eqs. ([H ]) -(fI5 ]) ]. The 
allowed electron-hopping processes out of the red-green 
bond dimers are for green electrons to hop along the c- 
axis bonds, becoming blue as they do so, and for red 
electrons to hop along the single &-axis bond, also be- 
coming blue. The change of color on hopping means that 




FIG. 7: (color online) Superexchange limit, (ac) plaquette, 
showing selected (a) horizontal and (b) "vertical" (ss/ot) 
dimer pairs (magenta ellipses), along with a schematic repre- 
sentation of allowed electron hopping processes. The orbital 
triplet states arc +— in both panels. 

there is no effect of the spin state in blocking possible 
processes. 

For finite elements in the t matrix, the final state must 
be a pair of spin singlets on the c-axis bonds with green- 
blue color triplet character, as represented in Fig. 7(b). 
Somewhat surprisingly, there exists precisely one second- 
order process which achieves this, namely a single ex- 
change of two initially red electrons on the &-axis bond 
when both the other electrons on the cluster are green 
(this situation is represented in Fig. 7). There is thus 
no need to consider correlated hopping processes occur- 
ring on both c-axis bonds simultaneously (which are in 
fact available only for the T z = states of all four bonds 
on the plaquette), as these are of order t^/U 3 . Thus at 
second order the plaquctte-flipping matrix has the max- 
imally sparse and inhomogcneous form 

*+-,+- = \J-> 
tfiv.pcr = otherwise. (27) 

This rather extreme form is due not only to the strict 
constraints set by the available color combinations, but 
also to the strong breaking of translational symmetry en- 
coded in the choice of rhombic plaquette. 

Under these circumstances, it is necessary to consider 
the second-order contributions to the QDM v term not 
contained in the dimer energy itself. From fluctuation 
processes on the three interdimer bonds one obtains 





-v, 


v +Q.+a = 




vo+.o+ = 


— £ J 


"+-,+- = 


-J, 


^00, 00 = 


4 J > 


+ = 


_I / 

2 ' 


V0-,0- = 


— — / 

4 J > 


V-0,-0 = 


8 J ' 


V . = 


4 J > 



v^iv.pa = otherwise. 
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(a) (b) 



(c) £J = a(&+ /^) 

FIG. 8: (color online) Two plaquette coverings in the "colum- 
nar dimer" class which maximize the number of fiippable pla- 
quettes for (ss/ot) dimers on the triangular lattice. Shown 
are (a) a highly regular configuration with active plaquettes 
of all three orientations and (b) a more generic configuration 
of active plaquettes with two orientations plus a single "de- 
fect" of the third orientation. The color scheme represents 
the fact (c) that the ground state of each active plaquette is 
the symmetrized combination of states involving dimers with 
only two out of the possible three colors in the relevant orbital 
triplet; however, this triplet changes with the orientation of 
both the plaquette and the occupied bond. 



The color constraint on the hopping possibilities is there- 
fore sufficiently strong that the result is a diagonal v ma- 
trix. The non-uniformity in the diagonal values is due to 
the effect of the single b bond, from which one color con- 
figuration can profit maximally and one cannot profit at 
all, while the other combinations may use this to differing 
degrees. 

From the manifestly strong effects on the t and v ma- 
trices arising from the choice of four-site plaquette units, 
one might conclude that symmetry-restoration effects 
due to summing over all of the available plaquettes and 
orientations are of prime importance, and that in this 
sense the consideration of a QDM does not advance the 
general understanding beyond that obtained in Ref. [l[. 
In the multicolored QDM, this strong breaking of sym- 
metry will result in a selection process favoring only those 
dimer colors maximizing the contributions from quantum 
fluctuations, and hence in a system of + and — dimers 
on each active plaquette. 

Whether it is possible to select the preferred dimer 
colors depends on the definition of the system. As in 
Sec. Ill, for models with no dimer number-mixing terms 
and fixed numbers of dimer of each color, the color sectors 
cannot mix and the system is constrained to maximize 
its energy with the available number of colors of each 
dimer. While the (os/st) case, in which the dimer color 
corresponds to the physical spin, is indeed constrained by 
this type of criterion (below), in the (ss/ot) case it is more 
straightforward to mix bond orbital sectors, as this can 
be done without changing the net electron orbital color 
balance, and therefore a selection of unequal numbers of 
T z states is not restricted. Thus the effective QDM for 
the (ss/ot) case is a type of sector-switching two-color 
model with four matrix elements in Vij and one in £y . 

Such a model, with |u| > \t\ and, crucially, v < 0, is 
known to be in the "columnar" VBC state, the general 



term applied to the set of dimer coverings which maxi- 
mize the number of fiippable plaquettes. Because of the 
constraint that — h plaquettes are inactive for 6-bond 
processes, the maximal number of active plaquettes ob- 
tainable in this "two-color" model is not the conventional 
1/6 but only 1/12. The v and t terms select automati- 
cally the optimal color state of the spin singlets, and the 
ground state of each active plaquette is the symmetrized 
combination of "horizontal" and "vertical" dimer pairs 
with energy — 2J — \v\ — t. On the triangular lattice of 
N sites, the total number of valid plaquette states of this 
type is a highly (but not extensively - the entropy is pro- 
portional to the perimeter of the system) degenerate set, 
of which two examples arc shown in Fig. 8. That the 
ground states are VBCs which break translational sym- 
metry does not mean that the system has no positional 
fluctuations: rather, this type of configuration maximizes 
the sole kinetic (t) term present, but this is not strong 
enough to melt the static (^-driven) plaquette order. In 
this case, the QDM result is quite unambiguous: the ra- 
tio deduced for \t/v\ is well in the columnar phase, and 
not close to a boundary. 



B. Superexchange limit, physical rj 



Here the ground state of the system is composed of 
(os/st) spin triplets. From the available electronic hop- 
ping processes, represented schematically in Fig. 9, it is 
clear that there is no second-order t term which can flip 
the orientation of a pair of (os/st) dimers. In contrast 
to the (ss/ot) case, the calculation of matrix elements 
thus requires the consideration of coherent electronic pro- 
cesses occurring at fourth order in the electronic Hamil- 
tonian. For a meaningful comparison between fourth 
order t and v terms in the relevant QDM, second-order 
contributions to v of the type computed in the previous 
subsection are treated as a renormalization of J. Let J 
denote 4t A /(U - 3J H ) 3 , then 

Vfiiy.pa = otherwise. 

The spin-conservation symmetry (the absence of a mag- 
netic field is assumed) makes the v matrix not only uni- 
form but, in combination with the spin-pairing possibil- 
ities and the fact that electron hopping involves an au- 
tomatic orbital color change, completely diagonal. The 
t matrix has more variety, because some overlap is per- 
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FIG. 9: (color online) Superexchange limit, (ac) plaquette, 
showing selected (a) horizontal and (b) "vertical" (os/st) 
dimer pairs (turquoise ellipses), along with a schematic rep- 
resentation of allowed electron hopping processes. Pairs of 
S z = 1 triplets are depicted in both panels. 

mittcd within the sectors of fixed total spin, and thus 
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In this case, there is blocking of the possible electronic 
processes leading to dimer flips by the allowed spin states: 
only the all-up and all-down spin configurations can 
profit from every available process. That the element 

too, oo is not symmetrical with and t . is a 

consequence of the breaking of translational symmetry 
in the plaquette choice. 

In this case the system is not free to choose its color 
state, which would mean choosing its real magnetic state 
(in effect this defines another set of unmixable "topolog- 
ical" sectors characterized by the total S z of the dimer 
ensemble). Under the reasonable assumption that the 
ground state will be net nonmagnetic, the system will 
have equal numbers of + and — triplets, but the ratio 
N± : Nq may lie anywhere between 0:1 and 1:0. The fact 
that matrix elements exist for many plaquette processes 
which mix the positions of and +/— triplets suggests 
that it is most favorable for the system to preserve the 
option of fluctuating among all the possible states. If 
the system were to freeze, locally or globally, to either 
limiting case, the possibilities for dynamical fluctuations 
would be strongly reduced. 

Thus one has deduced a true, three-color QDM, albeit 
one with predominantly diagonal elements in the 9x9 t 
and v matrices and with a ratio t^^/^^ = — j for 
most components. This is again a system which is in a 
robust columnar phase. A columnar phase in a system 
with three dimer colors, all (for the sake of illustration) 




FIG. 10: (color online) A columnar dimer covering maximiz- 
ing the number of flippable plaquettes on the triangular lattice 
for (os/st) dimers. 

equally occupied, would be one with a one-dimensional 
nature (hard but flexible strings of dimers) of the type 
shown in Fig. 10. 

Both in this case and for (ss/ot) dimers, individual 
matrix elements of v are rather larger than those in t. 
This result reflects the property of the model that all 
fluctuations away from and back to the starting configu- 
ration may contribute to v, while only quite specific pairs 
of color-switching processes are able to contribute to a t 
term. This places the system in the regime |d| > \t\ which 
is generically unsuitable for resonant states, and hence is 
the dominant physics deciding the (static or dynamic) na- 
ture of the ground state; the type of static phase is in turn 
determined by the sign of v. While it is again necessary 
to consider the restoration of translational symmetries 
broken rather crudely at the four-site plaquette level, it 
is difficult to conceive of circumstances under which this 
would bridge the gap, which includes changing the sign 
of v, from the columnar phase to an RVB phase. 

To conclude the analysis of this and the previous sub- 
section, from effective QDMs it appears very likely that 
the tig superexchange models on the triangular lattice fa- 
vor "columnar" VBC states. In the (ss/ot) case this is a 
type of two-color model with plaquette-ordered ground 
states, while in the (os/st) case it is a three-color model 
with dimer-ordered ground states. These states break 
the symmetry both in real space (lattice translation and 
rotation) and in the space of dimer color. This result is a 
consequence of the fact that virtual fluctuations promot- 
ing dimer resonance are simply too restrictive (in "color" , 
be this the orbital type or the real spin) to be numerous 
enough to compete qualitatively with virtual fluctuations 
profiting from static order. 



C. Direct— exchange limit, r\ = 

As demonstrated in Sec. II, in this model (fTT)) the 
dimer overlap factor is a = 0. Thus the model is natu- 
rally color-conserving. The conventional QDM plaquette 
terms are simply t = = v and the appropriate model is 
obtained from 6-bond, 3-dimer loops of the type shown 
in Fig. 3 (for which there are two orientations). These 
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FIG. 11: (color online) Representation of the six-bond tri- 
angular loop fluctuation which is the minimal kinetic process 
permitted in the direct-exchange limit. 

contribute at sixth order in the electronic Hamiltonian, 
both to a v' term for processes returning to the starting 
configuration and to a t 1 term for correlated hopping pro- 
cesses in the same direction around the triangle, which 
lead to a net shift of the dimer positions. For clarity this 
is illustrated again in Fig. 11, with the dimcrs retain- 
ing their color code (which here, in contrast to the pre- 
ceding subsections, does denote their orbital state). By 
analogy to the procedure mentioned above, fourth-order 
contributions to maintaining the static dimer configura- 
tion can also be renormalized into the energy of the dimer 
coverings (whose extensive degeneracy is not affected by 
these) . By writing out the wave functions for the triangle 
states, 

\4) = ^(|l&t2H)-|lH2&t>)(|3fft4<7+) (31) 

-|3 fl 4,4fft))(|5rt6r;)-|5r4.6rt)), 
= ^(|26t36|)-|2H36t»(|4.9t5.g|) (32) 

-|4 5 | 5 ff t»(|6rt Irl) - |6r| lrt)), 

and proceeding in a manner analogous to the previous 
two subsections, one finds 

t/ = -576^, t> = ^. (33) 

As above, the difference in magnitude between v' and t' 
arises due to the prevalence of connected processes not 
changing the dimer configuration over those leading to 
the specific resonance loop in question. 

Thus once again this is a model strongly in the colum- 
nar limit: based on a predominant ordering pattern and 
profiting from limited fluctuations (occurring only on se- 
lected generalized plaquette units). The meaning of the 
term "columnar" in standard QDM literature is a config- 
uration maximizing the number of flippable plaquettes, 
which for the t'—v' model means the number of up- and 
down-oriented triangles. Examples arc shown in Fig. 12. 
This columnar ground state has three possible color com- 
binations, specifically those where any two colors are used 
in two triangles each, plus a triangle chirality, and hence 
is six-fold degenerate. It does not possess equal numbers 
of dimers of each color (the ratio is 1:1:2). A ground state 
of a macroscopic system may be expected to have a do- 
main structure, and all states with equal dimer numbers 
will present defective versions of this phase. 




FIG. 12: (color online) "Columnar" dimer coverings of the 
triangular lattice maximizing the number of flippable six-site 
triangular units. Two of the six possibilities are illustrated: 
in panel (a) the red and green dimers (directions) each partic- 
ipate in two loops and the more numerous blue dimers in one; 
in panel (b) is shown the inequivalent state which still has 
blue as the unique direction but has the opposite chirality of 
dimer locations on the the up- and down-oriented triangular 
units. 

The overall energy gain of this state of the QDM rel- 
ative to the original manifold (which has eyB = — J' /&) 
is e = (v — t)(ng) per plaquette, which in the present 
case of six-bond loops is e 6 = (v' - t')/12 ~ -73tf/U 5 
per bond. However, because of the higher-order nature 
of this correction, it is not able to alter the conclusion^ 
that the one-dimensional solution, which has energy 
eiD = — ! In 2 = —0.23105 J' per bond, appears to lie 
lower than all dimer states. It is worth reiterating in this 
context that the present considerations are restricted to 
the case when 77 is identically zero; as noted in Sec. II, 
finite values of 77 lead to the stabilization of particular 
static VBC configurations which promote energy gains 
of 0(t7 3 ),— also insufficient to redress the energy bal- 
ance. The prevalence of the one-dimensional solution 
in a geometry of such high connectivity is testimony to 
the enormous degree of orbital-induced frustration in the 
full physical system. 

V. SUMMARY 

This manuscript pursues, by mapping to a set of quan- 
tum dimer nodels (QDMs), the nature of the ground 
state of the ti g spin-orbital model on the triangular lat- 
tice. The electronic Hamiltonian, investigated in detail in 
Ref. [if, represents an insulating 3d 1 electron system on 
the (111) planes of a system such as NaTiC>2, where the 
cubic structural symmetry of edge-sharing metal-oxygen 
octahedra leads to an unbroken, threefold orbital degen- 
eracy. This model was showni to have a very strong pref- 
erence for dimer-based states of no long-ranged magnetic 
or orbital order over its entire parameter range in both 
the ratios J'/ J (of the direct and superexchange contri- 
butions to the magnetic interactions) and Jh/U (of the 
Hund coupling to the on-site Coulomb repulsion). 

These dimer states were found to span a range of be- 
havior from highly resonant in the superexchange limit 
to quasi-static in the direct-exchange limit. Both cases 
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present fundamental challenges in determining the na- 
ture of the true ground state: the former is a candidate 
resonating valence-bond (RVB) state and the latter a sit- 
uation where subtle selection effects, sometimes known 
as "order-by-disorder," are responsible for choosing the 
true ground state from a very highly degenerate mani- 
fold. In Ref. Q only energetic studies were performed, 
which were unable to answer the topological questions as- 
sociated with the formation of an RVB state or to resolve 
the possible differences among quasi-static valence-bond 
coverings. 

These questions are addressed by mapping the elec- 
tronic Hamiltonian to the minimal QDM for each of three 
cases: (a) the pure superexchange limit for low Jh/U, 
(b) the pure superexchange limit for intermediate Jh/U, 
and (c) the direct-exchange limit for Jh/U — 0. In 
case (a), the ground state is composed of spin-singlet, 
orbital-triplet (ss/ot) dimer entities and the effective 
QDM has three dimer "colors" corresponding to their 
orbital state. In case (b), the ground state is composed 
of orbital-singlet, spin-triplet (os/st) dimers and the ef- 
fective QDM has again three colors, which correspond to 
the triplet spin components of the dimers. In case (c), 
there are three dimer colors corresponding to the orbital 
colors active in each of the three bond directions of the 
triangular lattice, but because these are locked to each 
other, the color is not a degree of freedom and the effec- 
tive QDM has one color. The "non-orthogonality catas- 
trophe" which affects conventional spin dimers, namely 
that all dimer coverings have finite overlap, is shown to 
be very strongly reduced or completely eliminated in the 
presence of a dimer color in such models. 

In order to analyze the properties of these specific 
QDMs, first the general multicolored QDM must be un- 
derstood. This is investigated for the two-color case, 
with some additional consideration of three-color QDMs 
where these contain further physics. Both the potential 
term v for adjacent dimers and the kinetic term t for 
these to flip direction on the four-site plaquette they de- 
fine become n 2 x n 2 matrices, where n is the number of 
colors. With such a choice of matrix elements, there re- 
mains only one Rokhsar-Kivelson (RK) point, where the 
equally weighted superposition of all possible dimer cov- 
erings (in a given topological sector) is an eigenstate, but 
the regime of parameters in which this (gapped, RVB) 
state is the ground state is expanded considerably. Non- 
trivial topological sectors are found for specific choices of 
matrix elements, notably those conserving the total num- 
ber of dimers of each color and those which allow this to 
be altered not in single steps but in units of two. This 
gives rise to topological excitations related to the color 
degree of freedom, which arc termed "color visons." 

The matrix elements and v^ UtPa are deduced 

from the electronic Hamiltonian for cases (a) and (b). 
The restrictions on hopping of electrons of different or- 
bital color and the breaking of translational symmetry 
contained in the consideration of four-site plaquettes lead 



in case (a) to very sparse and asymmetric matrices. This 
asymmetry selects only the T z = ±1 orbital triplet states, 
and the model is reduced to one with two colors. In case 
(b), the physical spin defines "topological" sectors which 
cannot mix. In both cases, elements of the v matrix are 
generally rather larger than those of the t matrix, and 
have in addition a negative sign. Thus the mapping to 
a QDM indicates strongly that the ground state of both 
models is a "columnar" plaquette phase, one which is 
based on a small subset of static dimer coverings which 
maximize the number of flippable four-site plaquettes: 
this number is 1/12 of the total in case (a) and 1 /6 in case 
(b). While these states gain energy from virtual dimer 
flipping processes (quantum fluctuations), such processes 
arc not sufficiently strong that they can "melt" the pla- 
quette or dimer order in favor of a completely resonant, 
symmetry-restored phase. 

In case (c), the overlap matrix elements, and t and 
v elements, on a four-site plaquette are identically zero. 
The leading quantum fluctuation processes for dimer res- 
onance occur on six-site, triangular units containing one 
dimer of each bond direction, and hence the correspond- 
ing minimal QDM is a t'-v' model defined on these "pla- 
quettes" . Once again, v' is negative and larger in mag- 
nitude than t', suggesting that the ground state of the 
system is again a generalized columnar covering (mean- 
ing one which maximizes the number of flippable trian- 
gles). These are illustrated and shown to have a six-fold 
degeneracy, corresponding to the bond direction and the 
triangle chirality. 

In summary, effective QDMs are used to obtain indi- 
cations as to the nature of the ground state of a complex 
electronic Hamiltonian for a spin-orbital model which is 
known to have no conventional magnetic or orbital order. 
In all cases of most interest, the ground state is found 
to be a resonance-stabilized spin-orbital VBC covering 
which maximizes quantum fluctuation processes on a re- 
stricted set of plaquettes, but does not allow the melting 
of a preferred order by these fluctuations. While the 
derivation of a QDM involves rather crude symmetry- 
breaking at the level of four- (or six-)site plaquettes, and 
the symmetry-restoring effects of fluctuations across all 
available plaquettes must be borne in mind, the columnar 
phases deduced from all three QDMs are quite robust, in- 
dicating that there is no nearby liquid phase. 
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